Tar mat formation prediction in late-charge reservoirs

ABSTRACT

A downhole tool, surface equipment, and/or remote equipment are utilized to obtain data associated with a subterranean hydrocarbon reservoir, fluid contained therein, and/or fluid obtained therefrom. At least one condition indicating that a density inversion exists in the fluid contained in the reservoir is identified from the data. Molecular sizes of fluid components contained within the reservoir are estimated from the data. A model of the density inversion is generated based on the data and molecular sizes. The density inversion model is utilized to estimate the density inversion amount and depth and time elapsed since the density inversion began to form within the reservoir. A model of a gravity-induced current of the density inversion is generated based on the data and the density inversion amount, depth, and elapsed time.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation-in-part of U.S. application Ser. No. 14/927,305, entitled “Tar Mat Formation Prediction in Late-Charge Reservoirs,” filed Oct. 29, 2015, which claims priority to and the benefit of U.S. Provisional Application No. 62/028,213, entitled “Prediction of Tar Mat Formation in Reservoirs with Late Light Hydrocarbons Charge,” filed Jul. 23, 2014. This application also claims priority to and the benefit of U.S. Provisional Application No. 62/234,433, entitled “Tar Mat Formation Prediction in Late-Charge Reservoirs,” filed Sep. 29, 2015. The entire disclosure of each of said applications is hereby incorporated herein by reference.

BACKGROUND OF THE DISCLOSURE

Tar mat formation studies aid understanding of the formation mechanism, locations, and intervals favorable for tar formation, which may be utilized to minimize drilling in tar-rich zones and/or avoid completing a tar zone. Thus, operators may take measures to understand the origin and distribution of tar mat in reservoirs, perhaps including providing a reliable prediction of specific occurrences of tar mat prior to drilling.

SUMMARY OF THE DISCLOSURE

This summary is provided to introduce a selection of concepts that are further described below in the detailed description. This summary is not intended to identify indispensable features of the claimed subject matter, nor is it intended for use as an aid in limiting the scope of the claimed subject matter.

The present disclosure introduces a system that includes a downhole tool operable within a wellbore extending from a wellsite surface into a subterranean hydrocarbon reservoir. The downhole tool is operable to obtain fluid from the reservoir and data associated with the obtained fluid. The system also includes an apparatus operable to identify at least one condition indicating that a density inversion exists in the fluid contained in the reservoir based on the obtained data.

The present disclosure also introduces a system that includes a downhole tool operable to obtain fluid from a subterranean reservoir at multiple depths within the reservoir. The downhole tool is also operable to obtain data associated with the obtained fluid, including gas-oil ratio (GOR) and composition data related to the fluid obtained at each of the depths. The system also includes an apparatus operable to identify at least one condition indicating that a density inversion exists in the fluid contained in the reservoir based on the obtained data. The at least one condition includes at least one of: (i) a first condition in which there exists a disequilibrium gradient of GOR relative to depth within the reservoir, as determined based on the GOR data at each of the plurality of depths; and/or (ii) a second condition in which there exists an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths.

These and additional aspects of the present disclosure are set forth in the description that follows, and/or may be learned by a person having ordinary skill in the art by reading the materials herein and/or practicing the principles described herein. At least some aspects of the present disclosure may be achieved via means recited in the attached claims.

BRIEF DESCRIPTION OF THE DRAWINGS

The present disclosure is best understood from the following detailed description when read with the accompanying figures. It is emphasized that, in accordance with the standard practice in the industry, various features are not drawn to scale. In fact, the dimensions of the various features may be arbitrarily increased or reduced for clarity of discussion.

FIG. 1 is a chart demonstrating one or more aspects related to the present disclosure.

FIG. 2 is a graph demonstrating one or more aspects related to the present disclosure.

FIG. 3 is a graph demonstrating one or more aspects related to the present disclosure.

FIG. 4 is a flow-chart diagram of at least a portion of a workflow and/or other method according to one or more aspects of the present disclosure.

FIG. 5 is a schematic view of at least a portion of apparatus according to one or more aspects of the present disclosure.

FIG. 6 is a schematic view of at least a portion of apparatus according to one or more aspects of the present disclosure.

FIG. 7 is a schematic view of at least a portion of apparatus according to one or more aspects of the present disclosure.

FIG. 8 is a schematic view of at least a portion of apparatus according to one or more aspects of the present disclosure.

FIG. 9 is a schematic view of at least a portion of apparatus according to one or more aspects of the present disclosure.

FIG. 10 is a flow-chart diagram of at least a portion of a method according to one or more aspects of the present disclosure.

FIG. 11 is a flow-chart diagram of at least a portion of a method according to one or more aspects of the present disclosure.

DETAILED DESCRIPTION

It is to be understood that the following disclosure provides many different embodiments, or examples, for implementing different features of various embodiments. Specific examples of components and arrangements are described below to simplify the present disclosure. These are, of course, merely examples and are not intended to be limiting. In addition, the present disclosure may repeat reference numerals and/or letters in the various examples. This repetition is for simplicity and clarity, and does not in itself dictate a relationship between the various embodiments and/or configurations discussed. Moreover, the formation of a first feature over or on a second feature in the description that follows may include embodiments in which the first and second features are formed in direct contact, and may also include embodiments in which additional features may be formed interposing the first and second features, such that the first and second features may not be in direct contact.

The present disclosure introduces workflows, methods, and/or other processes in which reservoir testing data is performed to determine parameters pertaining to a subterranean hydrocarbon reservoir. An asphaltene concentration gradient analysis (ACGA) workflow may also be performed to obtain information pertaining to molecular sizes of hydrocarbon components located in the reservoir. A density inversion model may be generated utilizing the determined parameters and the obtained molecular size information. The density inversion model may be utilized to determine an amount of density inversion and its appearance location and time. A gravity current model may be generated utilizing the density inversion model information. Downslope velocity and distance of denser fluid may then be determined based on the gravity current model, from which tar mat formation may be predicted. The present disclosure also introduces systems, apparatus, and program products for performing at least a portion of such workflows, methods, and/or other processes.

With respect to fields having a long and/or thick tar mat, it has been found that biodegradation, water-washing, and other mechanisms may not be able to provide a sufficient explanation and/or characterization of the tar mat. The present disclosure introduces aspects related to the tar mat formation mechanism due to late light hydrocarbons charge as an avenue to explaining the tar mat formation.

For example, the tar mat formation mechanism due to late light hydrocarbons (gas) charge in reservoirs can be utilized to qualitatively explain the tar mat phenomenon. Such an approach may account for the effects of late gas charge on fluid properties, such as gas-oil ratio (GOR) and asphaltene instability, which may demonstrate denser fluid occurrence due to the asphaltene cluster integration, and its migration to the reservoir base driven by gravity.

FIG. 1 is a schematic view depicting four stages, in geologic time (increasing from left to right relative to the page), of a reservoir 10 undergoing tar mat formation due to a late charge of light hydrocarbons 12. In the first depicted stage 20, before commencement of the charge 12 into the reservoir 10, the fluid in the reservoir 10 is assumed to be homogeneous fluid primarily including two main constituents, namely, a first constituent primarily comprising maltenes (perhaps comprising gas), and a second constituent primarily comprising asphaltenes.

In the second depicted stage 22 after the inflow of the hydrocarbon charge 12, the late (new) gas migrates rapidly to the top 14 of the reservoir 10, such as through channels (not shown) extending within the reservoir 10 in a direction having at least some vertical component. After migrating to the top of the reservoir 10 and creating a contact displacement 16, the gas dissolves into the oil phase under high pressure. Thus, the rich gas near the top 14 of the reservoir 10 creates a concentration difference 18 of gas in the oil column of the reservoir 10, extending from the contact displacement 16 down to a diffusion front 19.

Subsequently, as in the third depicted stage 24, the dissolved gas diffuses down along the oil column, which in turn causes a concentration difference of maltenes and asphaltenes. Thus, these three components diffuse down towards the base 15 of the reservoir 10. According to the solubility parameter of different components, the light gas components can expel the heavy asphaltenes components out the oil phase. Therefore, asphaltenes form larger clusters or otherwise become unstable and separate from the oil phase. These clusters aggregate and are pushed downward by the diffusion front 19. Thus, an inversion 13 of the bulk density will eventually exist, where fluid near the diffusion front 19 will be of higher density than that at lower depths, due to an accumulation of asphaltenes caused by repulsion of the asphaltenes below the diffusion front 19. After sufficient time, gravity will drive the denser fluid as a current down towards the base 15, as indicated in FIG. 1 by reference numeral 28. Eventually, as in the fourth depicted stage 26, the gravity-induced current of the asphaltenes will accumulate in the base 15, thus forming a high viscous tar mat 29 over a very large area.

However, the above mechanism and other existing knowledge merely explain tar mat formation qualitatively. Field case validation may provide a better understanding regarding how the density inversion 13 appeared as a result of the late gas charge 12, and how the denser fluid migrated to the base 15 and deposited there to form the tar mat 29 within geologic time. The present disclosure introduces one or more aspects that may be utilized to identify tar mat formation possibility in hydrocarbon reservoirs based on the above-described tar mat formation resulting from late light hydrocarbon charge. By integrating the diffusion mechanism with a gravity current model, the likelihood of tar mat occurring may be determined, and a particular status may be further formed, thus permitting analysis of tar mat formation within geologic time and reservoir connectivity.

FIG. 2 is a graph depicting an example prediction result of the fluid density in the oil column at different times. The graph includes a first curve 30 representing bulk fluid density relative to depth within the original reservoir before the late gas charge, depicting the assumption that the bulk fluid density of the oil column is substantially constant relative to depth. A second curve 31 represents the bulk fluid density relative to depth at about 10,000 years after commencement of the late gas charge, and depicts a density inversion 41 at about 15 meters (m) from the top of the reservoir. Another curve 32 represents the bulk fluid density relative to depth after about 50,000 years, and depicts a density inversion 42 at about 35 m from the top of the reservoir. The next three curves 33, 34, and 35 represent elapsed times of 100,000 years, 200.00 years, and 300,000 years, respectively, with corresponding density inversions 43, 44, and 45 existing at about 45 m, about 65 m, and about 85 m, respectively. The remaining three curves 36, 37, and 38 represent elapsed times of 500,000 years, 700,000 years, and 1,000,000 years, respectively, with corresponding density inversions 46, 47, and 48 existing at about 105 m, about 125 m, and about 155 m, respectively. Thus, as depicted in the example shown in FIG. 2, the density inversion caused denser fluid to flow downward towards the base, driven by gravity. One may also deduce from FIG. 2 that, for example, the accumulated, denser fluid probably formed a tar mat in the reservoir base.

The present disclosure introduces modeling the mechanism of tar mat formation due to late light hydrocarbons charge by the combination of a density inversion model and a gravity-induced current model. One or more aspects of the methods introduced herein may be utilized to dynamically predict the fluid density in the oil column and, in turn, obtain a better understanding of the possibility of tar mat occurrence and its formation status within geological time. This may include utilizing diffusion equations and/or related information to describe the molecular migration caused by the concentration difference. Such approach may be utilized in areas such as oil and gas, chemicals, fluid mechanics, and other example implementations in which denser fluid forms and migrates to the base of the reservoir, including implementations in which the density inversion and gravity induce horizontal fluid flow. For example, even for horizontal flow in a porous medium, the gravity current model may be successfully applied.

One or more aspects of the present disclosure may also be integrated with wireline reservoir testing (WRT) and ACGA. For example, WRT may provide information such as reservoir depth, pressure, temperature, fluid composition, density, and viscosity, among other parameters and/or characteristics, and the molecular sizes of the compositional components of fluid obtained via WRT may be obtained via ACGA.

An example workflow introduced herein may comprise three parts: data acquisition, the density inversion calculation, and the gravity-induced current simulation. Data from three aspects may be collected for the parameters/inputs to a simulator. Subsequently, the fluid density in the oil column over geologic time may be successfully predicted. The subsequently obtained understanding of tar mat formation may then be delivered after analyzing the denser fluid transport in the reservoir.

The density inversion model may be a one-dimensional (1-D), single-phase model. The reservoir fluid may be described as a single-liquid phase comprising three components: gas (typically methane and perhaps some other light hydrocarbons), maltenes, and asphaltenes. That is, as described above, once the molecular concentration difference appears, the rich gas dissolves in the top of the oil column and begins diffusing downward along the oil column, and then maltenes and asphaltenes also diffuse downward.

The mass conservation of gas, asphaltenes, and maltenes may be described as set forth below in Equations (1)-(3).

$\begin{matrix} {{\frac{\partial\left( c_{G} \right)}{\partial t} + \frac{\partial\left( {c_{G}u_{G}} \right)}{\partial z}} = 0} & (1) \\ {{\frac{\partial\left( c_{A} \right)}{\partial t} + \frac{\partial\left( {c_{A}u_{A}} \right)}{\partial z}} = 0} & (2) \\ {{\frac{\partial\left( c_{M} \right)}{\partial t} + \frac{\partial\left( {c_{M}u_{M}} \right)}{\partial z}} = 0} & (3) \end{matrix}$ where:

-   subscripts G, A, and M represent gas, asphaltenes, and maltenes,     respectively; -   c is molarity (molar concentration) of each component; -   u is migration velocity of each component; -   t is time; and -   z is the space coordinate originating from a lower boundary, such as     the oil-water-contact (OWC) at the base of the reservoir.

The initial condition can be expressed as set forth below in Equations (4)-(6). c _(G)=0.0,t=0  (4) c _(A) =c _(A0) ,t=0  (5) c _(M) =c _(M0) ,t=0  (6) where:

-   c_(A0) is the initial molarity of asphaltenes; and -   c_(M0) is the initial molarity of maltenes.

At the fixed bottom of the oil column, a non-permeable boundary condition may be applied to each component, as set forth below in Equations (7)-(9). u _(G)=0.0,z=0  (7) u _(A)=0.0,z=0  (8) u _(M)=0.0,z=0  (9)

The top of the oil column (i.e., the gas-oil-contact (GOC)) is a moving boundary, which is non-permeable for maltenes and asphaltenes, but gas is saturated with oil at the GOC. Thus, assuming the location of the GOC is Z(t), the component migration velocities may be expressed as set forth below in Equations (10)-(12). c _(G)(t,Z(t))=k _(G) c _(M) ,z=Z(t)  (10)

$\begin{matrix} {{{u_{A}\left( {t,{Z(t)}} \right)} = \frac{d\;{Z(t)}}{d\; t}},{z = {Z(t)}}} & (11) \\ {{{u_{M}\left( {t,{Z(t)}} \right)} = \frac{d\;{Z(t)}}{d\; t}},{z = {Z(t)}}} & (12) \end{matrix}$ where k_(G) is solubility value (equilibrium constant) of gas in maltenes at specified temperature and pressure.

Since no maltenes are lost from the system, Z(t) satisfies Equations (13) and (14) set forth below. Z(0)=L  (13) ∫₀ ^(Z(t)) c _(G)(t,Z(t))dz=c _(M0) L  (14) where L is the original height of the oil column.

To introduce constitutive diffusion equations to the equations above, the reference frame of diffusion flux may be defined. The reference frame may include mass average, molar average, volume average, and solvent fixed, for example. In some implementations, the solvent fixed frame of reference may provide a simpler form of boundary conditions. Thus, the diffusion flux J of each component may be expressed as set forth below in Equations (15)-(17). J _(G) =c _(G)(u _(G) −u _(M))  (15) J _(A) =c _(A)(u _(A) −u _(M))  (16) J _(M)=0  (17)

Accordingly, Equations (1)-(3) may be rewritten as set forth below in Equations (18)-(20).

$\begin{matrix} {\frac{\partial\left( c_{G} \right)}{\partial t} = {- \frac{\partial\left( {J_{G} + {c_{G}u_{M}}} \right)}{\partial z}}} & (18) \\ {\frac{\partial\left( c_{A} \right)}{\partial t} = {- \frac{\partial\left( {J_{A} + {c_{A}u_{M}}} \right)}{\partial z}}} & (19) \\ {\frac{\partial\left( c_{M} \right)}{\partial t} = {- \frac{\partial\left( {c_{M}u_{M}} \right)}{\partial z}}} & (20) \end{matrix}$

Similarly, the boundary conditions of Equations (7)-(12) may be rewritten as set forth below in Equations (21)-(26). J _(G)=0.0,z=0  (21) J _(A)=0.0,z=0  (22) J _(M)=0.0,z=0  (23) c _(G)(t,Z(t))=k _(G) c _(M) ,z=Z(t)  (24) J _(A)(t,Z(t))=0,z=Z(t)  (25)

$\begin{matrix} {{{u_{M}\left( {t,{Z(t)}} \right)} = \frac{d\;{Z(t)}}{d\; t}},{z = {Z(t)}}} & (26) \end{matrix}$

A space variable ξ may be expressed as set forth below in Equation (27). ξ=∫₀ ^(Z) c _(M)(t,z)dz  (27)

Equations (18)-(26) can be reduced to the format of a classic diffusion equation by transforming the variables that are functions of t and z to functions of t and ξ. Such transformation may utilize the relationships set forth below in Equations (28) and (29) for each function ƒ(t, z).

$\begin{matrix} {\left( \frac{\partial{f\left( {t,\xi} \right)}}{\partial\xi} \right)_{t} = {{\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{t}/\left( \frac{\partial{\xi\left( {t,z} \right)}}{\partial z} \right)_{t}} = {\frac{1}{c_{M}}\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{t}}}} & {(28)} \\ {\left( \frac{\partial{f\left( {t,\xi} \right)}}{\partial t} \right)_{\xi} = {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial t} \right)_{z} + {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{t}\left( \frac{\partial{z\left( {t,\xi} \right)}}{\partial t} \right)_{\xi}}}} & \\ {= {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{z} - {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{t}\frac{\left( \frac{\partial{\xi\left( {t,z} \right)}}{\partial t} \right)_{z}}{\left( \frac{\partial{\xi\left( {t,z} \right)}}{\partial z} \right)_{t}}}}} & \\ {= {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{z} - {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{t}\frac{\left( {\int_{0}^{z}{\frac{\partial c_{M}}{\partial t}\ d\; z}} \right)_{z}}{c_{M}\left( {t,z} \right)}}}} & \\ {= {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{z} + {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{t}\frac{\left( {\int_{0}^{z}{\frac{\partial\left( {c_{M}u_{M}} \right)}{\partial z}\ d\; z}} \right)_{z}}{c_{M}\left( {t,z} \right)}}}} & \\ {= {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial t} \right)_{z} + {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{t}u_{M}}}} & {{~~~~~}(29)} \end{matrix}$

This leads to Equations (30) and (31) set forth below.

$\begin{matrix} {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial z} \right)_{t} = {c_{M}\left( \frac{\partial{f\left( {t,\xi} \right)}}{\partial\xi} \right)}_{t}} & (30) \\ {\left( \frac{\partial{f\left( {t,z} \right)}}{\partial t} \right)_{z} = {\left( \frac{\partial{f\left( {t,\xi} \right)}}{\partial t} \right)_{\xi} + {c_{M}{u_{M}\left( \frac{\partial{f\left( {t,\xi} \right)}}{\partial\xi} \right)}_{t}}}} & (31) \end{matrix}$

Thus, Equations (18)-(20) may be expressed utilizing t and ξ by applying Equations (30) and (31), as set forth below in Equations (32)-(34).

$\begin{matrix} {\left( \frac{\partial c_{G}}{\partial t} \right)_{\xi} = {{{c_{M}{u_{M}\left( \frac{\partial c_{G}}{\partial\xi} \right)}_{t}} - {c_{M}\left( \frac{\partial}{\partial\xi} \right)}_{t}} = {{- {c_{M}\left( \frac{\partial J_{G}}{\partial\xi} \right)}_{t}} - {c_{M}{c_{G}\left( \frac{\partial u_{M}}{\partial\xi} \right)}_{t}}}}} & (32) \\ {\left( \frac{\partial c_{A}}{\partial t} \right)_{\xi} = {{{c_{M}{u_{M}\left( \frac{\partial c_{A}}{\partial\xi} \right)}_{t}} - {c_{M}\left( \frac{\partial}{\partial\xi} \right)}_{t}} = {{- {c_{M}\left( \frac{\partial J_{A}}{\partial\xi} \right)}_{t}} - {c_{M}{c_{A}\left( \frac{\partial u_{M}}{\partial\xi} \right)}_{t}}}}} & (33) \\ {\left( \frac{\partial c_{M}}{\partial t} \right)_{\xi} = {{{c_{M}{u_{M}\left( \frac{\partial c_{M}}{\partial\xi} \right)}_{t}} - {c_{M}\left( \frac{\partial}{\partial\xi} \right)}_{t}} = {- {c_{M}^{2}\left( \frac{\partial}{\partial\xi} \right)}_{t}}}} & (34) \end{matrix}$

Multiplying Equation (32) by 1/c_(M), multiplying Equation (34) by c_(G)/c_(M) ², subtracting the former from the latter, and introducing a new composition variable S_(G)=c_(G)/c_(M) results in Equation (35) set forth below.

$\begin{matrix} {\left( \frac{\partial S_{G}}{\partial t} \right)_{\xi} = {- \left( \frac{\partial J_{G}}{\partial\xi} \right)_{t}}} & (35) \end{matrix}$

Similarly utilizing Equations (33) and (34) and another new composition variable S_(A)=c_(A)/c_(M) results in Equation (36) set forth below.

$\begin{matrix} {\left( \frac{\partial S_{A}}{\partial t} \right)_{\xi} = {- \left( \frac{\partial J_{A}}{\partial\xi} \right)_{t}}} & (36) \end{matrix}$

The constitutive equations may now be derived such that J_(G) and J_(A) may be expressed by S_(G) and S_(A) with a linear relationship. Such derivation may start from the Maxwell Stefan theory, which links the chemical potential gradient to relative velocities between components, as set forth below in Equation (37).

$\begin{matrix} {{\frac{x_{i}}{RT}{\nabla_{z}\mu_{i}}} = {- {\sum_{j = 1}^{N}\frac{x_{i}{x_{j}\left( {u_{i} - u_{j}} \right)}}{D_{ij}}}}} & (37) \end{matrix}$ where:

-   x_(i) is mole fraction of component i; -   R is the universal gas constant; -   T is the WRT-measured temperature of the reservoir; -   ∇ is the gradient operator; -   μ_(i) is chemical potential of component i; -   N is the number of components i; and -   D is the Maxwell Stefan diffusion coefficient between components i     and j, where D_(ij)=D_(ji).

For systems considered within the scope of the present disclosure, two independent Maxwell Stefan equations for gas and asphaltenes can be written as set forth below in Equations (38)-(41).

$\begin{matrix} {\mspace{79mu}{{\frac{x_{G}}{RT}{\nabla_{z}\mu_{G}}} = {{- \frac{x_{G}{x_{A}\left( {u_{G} - u_{A}} \right)}}{D_{GA}}} - \frac{x_{G}{x_{M}\left( {u_{G} - u_{M}} \right)}}{D_{GM}}}}} & (38) \\ {\mspace{79mu}{{\frac{x_{A}}{RT}{\nabla_{z}\mu_{A}}} = {{- \frac{x_{A}{x_{G}\left( {u_{A} - u_{G}} \right)}}{D_{GA}}} - \frac{x_{A}{x_{M}\left( {u_{A} - u_{M}} \right)}}{D_{AM}}}}} & (39) \\ {{\nabla_{z}\mu_{G}} = {{{{- {RT}}\frac{{\frac{x_{A}}{x_{G}}J_{G}} - J_{A}}{c_{t}D_{GA}}} - {{RT}\frac{\frac{x_{M}}{x_{G}}J_{G}}{c_{t}D_{GM}}}} = {{- {{RTJ}_{G}\left( {\frac{\frac{x_{A}}{x_{G}}}{c_{t}D_{GA}} + \frac{\frac{x_{M}}{x_{G}}}{c_{t}D_{GM}}} \right)}} + {{RT}\frac{J_{A}}{c_{t}D_{GA}}}}}} & (40) \\ {{\nabla_{z}\mu_{A}} = {{{{- {RT}}\frac{{\frac{x_{G}}{x_{A}}J_{A}} - J_{G}}{c_{t}D_{GA}}} - {{RT}\frac{\frac{x_{M}}{x_{A}}J_{A}}{c_{t}D_{AM}}}} = {{{RT}\frac{J_{G}}{c_{t}D_{GA}}} - {{RTJ}_{A}\left( {\frac{\frac{x_{G}}{x_{A}}}{c_{t}D_{GA}} + \frac{\frac{x_{M}}{x_{A}}}{c_{t}D_{AM}}} \right)}}}} & (41) \end{matrix}$

With the definition of J_(i) in Equations (15)-(17) and the relationship of D_(ij)=D_(ji), Equation (42) set forth below may be obtained.

$\begin{matrix} \begin{matrix} {\begin{bmatrix} {\nabla_{z}µ_{G}} \\ {\nabla_{z}µ_{A}} \end{bmatrix} = {\frac{RT}{c_{t}}\begin{bmatrix} {{{- \frac{x_{A}}{x_{G}}}\frac{1}{D_{GA}}} - {\frac{x_{M}}{x_{G}}\frac{1}{D_{GM}}}} & \frac{1}{D_{GA}} \\ \frac{1}{D_{GA}} & {{{- \frac{x_{G}}{x_{A}}}\frac{1}{D_{GA}}} - {\frac{x_{M}}{x_{A}}\frac{1}{D_{AM}}}} \end{bmatrix}}} \\ {\begin{bmatrix} J_{G} \\ J_{A} \end{bmatrix}} \\ {= {\frac{RT}{c_{t}}\begin{bmatrix} {{{- \frac{S_{A}}{S_{G}}}\frac{1}{D_{GA}}} - {\frac{1}{S_{G}}\frac{1}{D_{GM}}}} & \frac{1}{D_{GA}} \\ \frac{1}{D_{GA}} & {{{- \frac{S_{G}}{S_{A}}}\frac{1}{D_{GA}}} - {\frac{1}{S_{A}}\frac{1}{D_{AM}}}} \end{bmatrix}}} \\ {\left. \begin{bmatrix} J_{G} \\ J_{A} \end{bmatrix}\Leftrightarrow\begin{bmatrix} J_{G} \\ J_{A} \end{bmatrix} \right.} \\ {= {\frac{c_{t}}{RT}\begin{bmatrix} {{{- \frac{S_{A}}{S_{G}}}\frac{1}{D_{GA}}} - {\frac{1}{S_{G}}\frac{1}{D_{GM}}}} & \frac{1}{D_{GA}} \\ \frac{1}{D_{GA}} & {{{- \frac{S_{G}}{S_{A}}}\frac{1}{D_{GA}}} - {\frac{1}{S_{A}}\frac{1}{D_{AM}}}} \end{bmatrix}}^{- 1}} \\ {C_{M}\begin{bmatrix} {\nabla_{\xi}\mu_{G}} \\ {\nabla_{\xi}\mu_{A}} \end{bmatrix}} \\ {= {c_{M}^{2}\left( {1 + S_{A} + S_{G}} \right)}} \\ {\begin{bmatrix} {{{- \frac{S_{A}}{S_{G}}}\frac{1}{D_{GA}}} - {\frac{1}{S_{G}}\frac{1}{D_{GM}}}} & \frac{1}{D_{GA}} \\ \frac{1}{D_{GA}} & {{{- \frac{S_{G}}{S_{A}}}\frac{1}{D_{GA}}} - {\frac{1}{S_{A}}\frac{1}{D_{AM}}}} \end{bmatrix}^{- 1}} \\ {\frac{1}{RT}\begin{bmatrix} {\nabla_{\xi}\mu_{G}} \\ {\nabla_{\xi}\mu_{A}} \end{bmatrix}} \\ {= {{- c_{M}^{2}}\frac{1 + S_{A} + S_{G}}{{D_{AM}S_{G}} + D_{GA} + {D_{GM}S_{A}}}}} \\ {\begin{bmatrix} {D_{GM}{S_{G}\left( {{D_{AM}S_{G}} + D_{GA}} \right)}} & {D_{AM}D_{GM}S_{A}S_{G}} \\ {D_{AM}D_{GM}S_{A}S_{G}} & {D_{AM}{S_{A}\left( {D_{GA} + {D_{GM}S_{A}}} \right)}} \end{bmatrix}} \\ {\frac{1}{RT}\begin{bmatrix} {\nabla_{\xi}\mu_{G}} \\ {\nabla_{\xi}\mu_{A}} \end{bmatrix}} \end{matrix} & (42) \end{matrix}$

The flux equations may then be obtained, as set forth below in Equation (43).

$\begin{matrix} {\begin{bmatrix} J_{G} \\ J_{A} \end{bmatrix} = {{\frac{- {c_{M}^{2}\left( {1 + S_{A} + S_{G}} \right)}}{{D_{AM}S_{G}} + D_{GA} + {D_{GM}S_{A}}}\left\lbrack \begin{matrix} {D_{GM}{S_{G}\left( {{D_{AM}S_{G}} + D_{GA}} \right)}} & {D_{AM}D_{GM}S_{A}S_{G}} \\ {D_{AM}D_{GM}S_{A}S_{G}} & {D_{AM}{S_{A}\left( {D_{GA} + {D_{GM}S_{A}}} \right)}} \end{matrix} \right\rbrack}{{\frac{1}{RT}\begin{bmatrix} \frac{\partial\mu_{G}}{\partial S_{G}} & \frac{\partial\mu_{G}}{\partial S_{A}} \\ \frac{\partial\mu_{A}}{\partial S_{G}} & \frac{\partial\mu_{A}}{\partial S_{A}} \end{bmatrix}}\begin{bmatrix} {\nabla_{\xi}S_{G}} \\ {\nabla_{\xi}S_{A}} \end{bmatrix}}}} & (43) \end{matrix}$

Based on regular solution theory, the chemical potential of each component i is given by Equation (44) set forth below.

$\begin{matrix} {\mu_{i} = {\mu_{i}^{0} + {{RT}\left\{ {{\ln\frac{x_{i}\nabla_{i}}{V}} + 1 - \frac{\nabla_{i}}{V} + {\frac{\nabla_{i}}{RT}{\sum_{j}{\sum_{k}{\phi_{j}{\phi_{k}\left\lbrack {\left( {\delta_{i} - \delta_{j}} \right)^{2} + {2\; l_{ij}\delta_{i}\delta_{j}} - {0.5\left( {\delta_{j} - \delta_{k}} \right)^{2}} - {l_{jk}\delta_{j}\delta_{k}}} \right\rbrack}}}}}} \right\}}}} & (44) \end{matrix}$ where:

-   V _(i) is the partial molar volume of component i; -   φ is the volume fraction of component i; -   δ_(i) is the solubility parameter of component i; -   l_(ij) is the binary interaction parameter between components i and     j; and -   an N component mixture solubility parameter δ can be determined     utilizing Equations (45) and (46) set forth below.

$\begin{matrix} {\delta = {{\sum_{i}{\phi_{i}\delta_{i}}} = \frac{\sum_{i}{x_{i}{\overset{\_}{v}}_{i}\delta_{i}}}{\sum_{j}{x_{j}{\overset{\_}{v}}_{j}}}}} & (45) \end{matrix}$ V=Σ _(i) x _(i) V _(i)  (46)

where:

-   V is the molar volume of the mixture.

The chemical potential in Equation (44) can also be differentiated with respect to x_(i). Thus, Equations (47) and (48) set forth below may be obtained.

$\begin{matrix} {\mspace{79mu}{\frac{\partial x_{i}}{\partial S_{j}} = {\frac{\partial\left( \frac{S_{i}}{1 + {\sum_{k = 1}^{N - 1}S_{k}}} \right)}{\partial S_{j}} = \frac{{\left( {1 + {\sum_{k = 1}^{N - 1}S_{k}}} \right)\Delta_{ij}} - S_{i}}{\left( {1 + {\sum_{k = 1}^{N - 1}S_{k}}} \right)^{2}}}}} & (47) \\ {\frac{\partial\phi_{i}}{\partial S_{j}} = {\frac{\partial\left( \frac{{\overset{\_}{V}}_{i}S_{i}}{{\overset{\_}{V}}_{N} + {\sum_{k = 1}^{N - 1}{{\overset{\_}{V}}_{k}S_{k}}}} \right)}{\partial S_{j}} = \frac{{\left( {{\overset{\_}{V}}_{N} + {\sum_{k = 1}^{N - 1}{{\overset{\_}{V}}_{k}S_{k}}}} \right){\overset{\_}{V}}_{j}\Delta_{ij}} - S_{i}}{\left( {{\overset{\_}{V}}_{N} + {\sum_{k = 1}^{N - 1}{{\overset{\_}{V}}_{k}S_{k}}}} \right)^{2}}}} & (48) \end{matrix}$ where:

-   Δ_(ij) is the Kronecker delta.

The models described above can be further simplified by assuming that cross-term diffusion effects can be neglected, as set forth below in Equation (49).

$\begin{matrix} {\begin{pmatrix} \frac{\partial x_{G}}{\partial t} \\ \frac{\partial x_{A}}{\partial t} \end{pmatrix} = {\nabla\left\lbrack {\begin{pmatrix} {x_{G}D_{G}^{eff}} & 0 \\ 0 & {x_{A}D_{A}^{eff}} \end{pmatrix}\frac{1}{RT}\begin{pmatrix} \frac{\partial\mu_{G}}{\partial x_{G}} & \frac{\partial\mu_{G}}{\partial x_{A}} \\ \frac{\partial\mu_{A}}{\partial x_{G}} & \frac{\partial\mu_{A}}{\partial x_{A}} \end{pmatrix}\begin{pmatrix} \frac{\partial x_{G}}{\partial z} \\ \frac{\partial x_{A}}{\partial z} \end{pmatrix}} \right.}} & (49) \end{matrix}$ The initial conditions may be expressed as set forth below in Equation (50). x _(i)(0,z)=x _(i) ⁰ ,i=G,A  (50) The boundary conditions at the bottom of the hydrocarbon column can be expressed as set forth below in Equation (51).

$\begin{matrix} {{\frac{\partial{x_{i}\left( {t,0} \right)}}{\partial z} = 0},{i = G},A} & (51) \end{matrix}$ The boundary conditions at the GOC can be expressed as set forth below in Equation (52).

$\begin{matrix} {{\frac{\partial{x_{i}\left( {t,{z(t)}} \right)}}{\partial z} = 0},{i = G},A} & (52) \end{matrix}$ The location of the GOC can be approximated as set forth below in Equation (53).

$\begin{matrix} {{Z(t)} = {{Z(0)} + {\int_{0}^{t}{\frac{J_{G}^{GOC}}{C_{G}^{{GOC} - O} - C_{G}^{{GOC} - G}}\ d\; t}}}} & (53) \end{matrix}$ where:

-   J_(G) ^(GOC) is the gas mass flux charging into the oil column; -   C_(G) ^(GOC-O) is the molar concentration of gas at the GOC on the     oil side; and -   C_(G) ^(GOC-G) is the molar concentration of gas at the GOC on the     gas side.

Accordingly, the density inversion may be determined utilizing Equations (49)-(53), obtained via ones of Equations (1)-(48) or solving ones of Equations (1)-(48) directly. Table 1 lists examples of the inputs. Various ones of the parameters may be estimated (via WRT or otherwise) prior to the actual determination.

TABLE 1 Example Density Inversion Model Inputs and Origins Parameter Example Origin Reservoir Pressure P_(res) WRT Reservoir Temperature T_(res) WRT Reservoir Fluid Density ρ WRT Reservoir Fluid Viscosity μ WRT Reservoir Fluid Composition WRT Reservoir Depth WRT Time Setting Reservoir Fluid Molecular Size ACGA D_(GA) Estimate D_(GM) Estimate D_(AM) Estimate

As described above, the molecular sizes of the fluid components may be estimated via ACGA utilizing data obtained via WRT, such as: density, optical density, a bulk fluid solubility parameter, and molar volume of the fluid obtained via WRT at multiple depths within the reservoir; molar volume, partial density, and a compositional component solubility parameter of at least one of the compositional components of the WRT-obtained fluid at the various depths; and temperature of the reservoir. For example, estimating the molecular sizes utilizing ACGA and the WRT-obtained data may be via Equation (54) set forth below.

$\begin{matrix} {{\frac{{OD}\left( h_{2} \right)}{{OD}\left( h_{1} \right)} = {\exp\left\{ {\frac{v_{a}{g\left( {\rho - \rho_{a}} \right)}\left( {h_{2} - h_{1}} \right)}{RT} + {\frac{v_{a}}{RT}\left\lbrack {\left( {\delta_{a} - \delta} \right)_{h_{1}}^{2} - \left( {\delta_{a} - \delta} \right)_{h_{2}}^{2}} \right\rbrack} + \left\lbrack {\left( \frac{v_{a}}{v} \right)_{h_{2}} - \left( \frac{v_{a}}{v} \right)_{h_{1}}} \right\rbrack} \right\}}},} & (54) \end{matrix}$ where:

-   OD(h₁) is measured optical density of the WRT-obtained fluid at a     first depth h₁; -   OD(h₂) is measured optical density of the WRT-obtained fluid at a     second depth h₂; -   ν_(a) is estimated molar volume of the asphaltene portion of the     fluid contained within the reservoir; -   ν is estimated molar volume of the fluid contained within the     reservoir; and -   g is gravitational acceleration.

From the density inversion model described above, the density inversion amount and its variation with depth and time may be obtained. Subsequently, tar mat formation may be predicted utilizing a gravity current model. The gravity current model may be a two-dimensional (2-D) model describing the gravity current propagation in an inclined porous layer, such as illustrated in the schematic view depicted in FIG. 3.

Combining Darcy's law, the conservation of mass, and the processes of convection and diffusion, the denser fluid can be described as set forth below in Equation (55).

$\begin{matrix} {\frac{\partial H}{\partial t} = {\frac{{\Delta\rho}\;{gk}}{\phi\mu}\left( {{{–sin\theta}\frac{\partial H}{\partial x}} + {\frac{\cos\;\theta}{2}{\nabla^{2}H^{2}}}} \right)}} & (55) \end{matrix}$ where:

-   H is height of the gravity current at given location; -   Δρ is the density inversion; -   k is permeability of the layer; -   φ is porosity; -   μ is fluid viscosity; -   x is distance from the initial influx location; and -   θ is dip angle of the slope.

The first term on the right-hand side of Equation (55) represents the gravity-driven convection downslope, and the second term represents spreading perpendicular to the slope. It is noted that the gravity current downslope velocity may be of more importance relative to velocities in other directions in tar mat studies, because the tar mat formation time is more dependent on the downslope slipping. By simulating the gravity current within geological time, the downslope velocity and the experienced distance of denser fluid in specific time can be determined or estimated.

FIG. 4 is a flow-chart diagram of at least a portion of an example implementation of a method (50) incorporating the aspects described above. The method (50) may comprise performing (52) WRT, such as may utilize one or more downhole tools having one or more aspects in common with the wireline downhole tools and/or apparatus described below. Performing (52) the WRT may also or instead utilize downhole tools and/or apparatus conveyed downhole via means other than wireline, such as drillstring, coiled tubing, TLC tubing, and/or other conveyance means. The data obtained via performing (52) the WRT and/or other means may be utilized to calculate, estimate, predict, determine, and/or otherwise obtain (53) depth, mass composition, P_(res), T_(res), ρ, μ, and/or other parameters and/or characteristics of the reservoir.

The method (50) also comprises estimating (55) the molecule sizes of the components of the fluid contained within the reservoir. As described above, estimating (55) the molecule sizes may utilize ACGA, such as via Equation (54) set forth above.

The obtained (53) reservoir/fluid parameters and the estimated (55) component molecule sizes may then be utilized to generate (56) the density inversion model, such as via utilization of Equations (49)-(53), obtained via ones of Equations (1)-(48) or ones of Equations (1)-(48) directly.

The generated (56) density inversion model may then be utilized to estimate (57) the density inversion amount, location, and elapsed time since forming (e.g., since the start of the late gas charge). This information may then be utilized to model (58) the resulting gravity-induced current, such as via Equation (55) set forth above. The downslope velocity and distance of denser fluid may then be estimated (59) utilizing the gravity-induced current model.

The example method (50) depicted in FIG. 4, and well as other implementations within the scope of the present disclosure, may be utilized to calculate the fluid density inversion and simulate its following transport in a porous medium, and may be applied to explain the real tar mat formation process in subterranean reservoirs. The simulation results, including the density inversion and its downslope velocity and distance, may subsequently be utilized to analyze how likely it is that tar mat has formed, as well as where and when tar mat occurs and its development within geological time. Such information may aid in understanding the tar mat formation history and future process. Additionally, the reservoir conditions that may be likely to result in tar mat can be recognized. Accordingly, drilling and production resources may be utilized more effectively, such as by reducing or avoiding the risk of drilling into and/or completing a tar mat zone.

It is also noted that modeling the density inversion within the scope of the present disclosure is not limited to the example described above. For example, the following description relates to another example implementation for modeling the density inversion within the scope of the present disclosure. One or more aspects of the example implementation described above may also be utilized in combination with one or more aspects of the example implementation described below, among other example implementations also within the scope of the present disclosure.

As with the example described above, the following model is a 1-D, three-component, single-phase model. The reservoir fluid is described as a single fluid phase limited to three components: gas (typically methane and perhaps some other light hydrocarbons); maltenes (live oil components excluding asphaltenes); and asphaltenes. Each component is again represented by G, M, and A, respectively. The original gas in the oil is not explicitly included in the analysis, and may be considered to be included with the maltenes. A preexisting gas cap at the top of the oil reservoir is also not part of this model.

According to the conservation of mass for each component in a unit volume of fluid, the rate of change of concentration in the volume equals the net flux into this volume. As the model is 1-D, the mass conservation equation for each component may be as set forth below in Equation (56).

$\begin{matrix} {\frac{\partial c_{i}}{\partial t} = {- \frac{\partial J_{d}}{\partial z}}} & (56) \end{matrix}$ The diffusion flux J_(d) is described by Fick's first law, as set forth below in Equation (57).

$\begin{matrix} {J_{d} = {{- D}\frac{\partial c_{i}}{\partial z}}} & (57) \end{matrix}$

For asphaltenes, in addition to the diffusion, reducing solubility of asphaltenes in high solution gas oil due to a gas charge is considered, as referred to above as asphaltene expulsion. A proposed numerical value of the solubility parameter δ as √{square root over ((ΔH−RT)/ν_(m))}. The product of the sum of volume fractions and the solubility parameter of each component can estimate the solubility parameter of mixtures. As solution gas concentration increases, the difference in solubility parameter between asphaltenes and the fluid (δ_(A)−δ_(fluid)) becomes large, leading to asphaltene expulsion from the fluid phase. In other words, the dissolving capability of the fluid for asphaltenes is weakened as more gas enters the system. The derived FHZ dynamic equation describes the dynamics of a multicomponent mixture under convection, molecular diffusion, solubility effect, and gravity effect. In this case, due to the gas charge, the solubility effect dominates. Therefore, only molecular diffusion and solubility are considered. Accordingly, the asphaltene expulsion flux J_(e,A) may be expressed as set forth below in Equation (58).

$\begin{matrix} {J_{e,A} = {{Kc}\frac{\partial c_{l}}{\partial x}}} & (58) \end{matrix}$

The parameter K represents the solubility effects of each component, such as the rate at which gas expels asphaltenes from their original location, and can also be expressed as set forth below in Equation (59).

$\begin{matrix} {K = {\alpha\frac{2V_{G}{V_{A}\left( {\delta_{M} - \delta_{G}} \right)}\left( {\delta_{A} - \delta_{M}} \right)}{RT}}} & (59) \end{matrix}$ where α is an adjustable constant that accounts for the effects of the compensation of the asphaltene solubility parameter and of stability from asphaltene diffusion. For example, α may be set as 1.

Combining Equations (56)-(59), the transport equations for the three component system in single-phase under a gas charge may be as set forth below in Equations (60)-(62).

$\begin{matrix} {\frac{\partial c_{G}}{\partial t} = {D_{G}\frac{\partial^{2}c_{G}}{\partial x^{2}}}} & (60) \\ {\frac{\partial c_{M}}{\partial t} = {D_{M}\frac{\partial^{2}c_{M}}{\partial x^{2}}}} & (61) \\ {\frac{\partial c_{A}}{\partial t} = {{D_{A}\frac{\partial^{2}c_{A}}{\partial x^{2}}} + {K\left( {{\frac{\partial c_{A}}{\partial x}\frac{\partial c_{G}}{\partial x}} + {c_{A}\frac{\partial^{\; 2}c_{G}}{\partial x^{2}}}} \right)}}} & (62) \end{matrix}$ where D is the effective diffusion coefficient independent of time and depth, as described above.

As above, the initial condition is defined as a homogenous mixture of maltenes and asphaltenes in the oil reservoir without charged gas, as set forth below in Equation (63). c _(G)(0,x)=0;c _(i)(0,x)=constant_(i) i=M,A  (63)

The charge of gas into the oil is considered to occur at the top of the oil reservoir. It is also assumed that the original homogenous mixture of maltenes and asphaltenes remains in the oil reservoir without leaving the system. Thus, the boundary conditions are as set forth below in Equations (64)-(66). c _(G)(t,top)=constant  (64)

$\begin{matrix} {{{\frac{\partial c_{i}}{\partial x}}_{x = {top}} = {{0\mspace{31mu} i} = M}},A} & (65) \\ {{{\frac{\partial c_{i}}{\partial x}}_{x = {base}} = {{0\mspace{31mu} i} = G}},M,A} & (66) \end{matrix}$

FIG. 5 is a schematic view of an example WRT system 100 that may be employed onshore and/or offshore according to one or more aspects of the present disclosure, representing an example environment in which one or more aspects described above may be implemented, such as to perform one or more aspects of the method (50) depicted in FIG. 4 and/or other methods and/or processes within the scope of the present disclosure. As depicted in FIG. 5, a downhole tool 105 may be suspended from a platform, rig, derrick, and/or other wellsite structure 110 in a wellbore 111 formed in one or more subterranean reservoirs 10. The downhole tool 105 may be or comprise one or more tools, one or more of which may be or comprise an acoustic tool, a conveyance tool, a density tool, an electromagnetic (EM) tool, a reservoir evaluation tool (also known in the art as a formation evaluation tool), a magnetic resonance tool, a monitoring tool, a neutron tool, a nuclear tool, a photoelectric factor tool, a porosity tool, a reservoir characterization tool, a resistivity tool, a sampling tool, a seismic tool, a surveying tool, and/or a telemetry tool, although other downhole tools are also within the scope of the present disclosure.

The downhole tool 105 may be deployed from the wellsite structure 110 into the wellbore 111 via a conveyance means 115, which may be or comprise a wireline cable, a slickline cable, and/or coiled tubing, although other means for conveying the downhole tool 105 within the wellbore 111 are also within the scope of the present disclosure. As the downhole tool 105 operates, outputs of various numbers and/or types from the downhole tool 105 and/or components thereof (one of which is designated by reference numeral 120) may be sent via the conveyance means 115 and/or otherwise to a logging and control system and/or other surface equipment 160 at surface, and/or may be stored in various numbers and/or types of memory for subsequent recall and/or processing after the downhole tool 105 is retrieved to surface. Such data may also be transmitted from the surface equipment 160 to remote equipment 101 disposed remote from the WRT system 100, perhaps hundreds or thousands of kilometers away. For example, such transmission may be via one or more networks 102, such as may include one or more of a cellular network, a satellite network, wide area networks, local area networks, and/or other types of networks, including wired and/or wireless networks.

FIG. 6 is a schematic view of an example implementation of the WRT system 100 shown in FIG. 5, including an example implementation of the downhole tool 105 that is designated in FIG. 6 by reference numeral 200. The downhole tool 200 is operable to engage a portion of a sidewall of the wellbore 111 penetrating the reservoir 10. The downhole tool 200 may be utilized to perform one or more aspects of the method (50) depicted in FIG. 4 and/or other methods and/or processes within the scope of the present disclosure.

The downhole tool 200 is suspended in the wellbore 111 from a lower end of the conveyance means 115. At the surface, the conveyance means 115 may be communicatively coupled to an electronics and processing system and/or other surface equipment 160. The surface equipment 160 may include a controller having an interface configured to receive commands from a surface operator. The surface equipment 160 may further comprise a processor configured to implement one or more aspects of the methods described herein.

The downhole tool 200 may comprise a telemetry module 210, a reservoir test module 214, and a sample module 226 (among other example implementations of the component 120 depicted in FIG. 5). Although the telemetry module 210 is shown as being implemented separate from the reservoir test module 214, the telemetry module 210 may be implemented in the reservoir test module 214. The downhole tool 200 may also comprise additional components at various locations, such as a module 208 above the telemetry module 210 and/or a module 228 below the sample module 226, which may have varying functionality within the scope of the present disclosure.

The reservoir test module 214 may comprise a selectively extendable probe assembly 216 and a selectively extendable anchoring member 218 that are respectively arranged on opposing sides of the downhole tool 200. The probe assembly 216 may be operable to selectively seal off or isolate selected portions of the sidewall of the wellbore 111. For example, the probe assembly 216 may comprise a sealing pad that may be urged against the sidewall of the wellbore 111 in a sealing manner to prevent movement of fluid into or out of the reservoir 10 other than through the probe assembly 216. The probe assembly 216 may thus be operable to fluidly couple a pump 221 and/or other components of the reservoir tester 214 to the adjacent reservoir 10. Accordingly, the reservoir tester 214 may be utilized to obtain fluid samples from the reservoir 10 by extracting fluid from the reservoir 10 using the pump 221. A fluid sample may thereafter be expelled through a port (not shown) into the wellbore 111, or the sample may be directed to one or more detachable chambers 227 disposed in the sample module 226. In turn, the detachable fluid collecting chambers 227 may receive and retain the reservoir fluid for subsequent testing at surface or a testing facility. The detachable sample chambers 227 may be certified for highway and/or other transportation. The module 208 and/or the module 228 may comprise additional sample chambers 227, which may also be detachable and/or certified for highway and/or other transportation.

The reservoir tester 214 may also be utilized to inject fluid into the reservoir 10 by, for example, pumping the fluid from one or more fluid collecting chambers disposed in the sample module 226 via the pump 221. Moreover, while the downhole tool 200 is depicted as comprising one pump 221, it may also comprise multiple pumps. The pump 221 and/or other pumps of the downhole tool 200 may also comprise a reversible pump operable to pump in two directions (e.g., into and out of the reservoir 10, into and out of the collecting chamber(s) of the sample module 226, etc.).

The probe assembly 216 may comprise one or more sensors 222 adjacent a port of the probe assembly 216, among other possible locations. The sensors 222 may be utilized to determine petrophysical parameters of a portion of the reservoir 10 proximate the probe assembly 216. For example, the sensors 222 may be configured to measure or detect one or more of pressure, temperature, composition, electric resistivity, dielectric constant, magnetic resonance relaxation time, nuclear radiation, and/or combinations thereof, although other types of sensors are also within the scope of the present disclosure.

The reservoir tester 214 may also comprise a fluid sensing unit 220 through which obtained fluid samples may flow, such as to measure properties and/or composition data of the sampled fluid. For example, the fluid sensing unit 220 may comprise one or more of a spectrometer, a fluorescence sensor, an optical fluid analyzer, a density and/or viscosity sensor, and/or a pressure and/or temperature sensor, among others.

The telemetry module 210 and/or another portion of the downhole tool 200 may comprise a downhole controller and/or control system 212 communicatively coupled to the surface equipment 160. The surface equipment 160 and/or the downhole controller and/or control system 212 may be configured to control the probe assembly 216 and/or the extraction of fluid samples from the reservoir 10, such as via the pumping rate of pump 221. The surface equipment 160 and/or the downhole controller and/or control system 212 may be further configured to analyze and/or process data obtained from sensors disposed in the fluid sensing unit 220 and/or the sensors 222, store measurements or processed data, and/or communicate measurements or processed data to surface or another component for subsequent analysis.

FIG. 7 depicts an example implementation of a spectrometer 305 and a control/monitoring system 310 that may form a portion of the downhole tool 100 shown in FIG. 5, the downhole tool 200 shown in FIG. 6, and/or other downhole tools within the scope of the present disclosure, that may be utilized to estimate or determine optical density (OD) and/or other properties of fluid obtained from the reservoir 10. The spectrometer 305 may comprise a light source 315 and a detector 320 disposed on opposite sides of a flowline 325 through which the pumped reservoir fluid flows, as indicated by arrow 330. The spectrometer 305 may be located at various possible locations along the flowline 325 that directs the reservoir fluid through the downhole tool. Although a single light source 315 is depicted in the example shown in FIG. 7, the spectrometer 305 may include additional light sources 315.

The detector 320 senses the light that passes through the reservoir fluid in the flowline 325. The detector 320 may include one or more detector elements 335 that may each be operable to measure the amount of light transmitted at a certain wavelength. For example, the detector elements 335 may detect the light transmitted from the visible to near-infrared within a range of 1, 5, 10, 20, or more different wavelengths ranging between about 400 nm and about 2200 nm. However, other numbers of wavelengths (corresponding to the number of detector elements) and other ranges of wavelengths are also within the scope of the present disclosure. For example, optical characteristics of the reservoir fluid may be detected at a range of wavelengths, such as the near infrared (NIR) wavelength range of approximately 800-2500 nm, 1500-2050 nm, or 1600-1800 nm. Estimations of reservoir fluid properties according to one or more aspects of the present disclosure may utilize optical data collected at a single wavelength, at multiple wavelengths, a range of wavelengths, and/or multiple ranges of wavelengths.

The spectrometer 305 may measure one or more optical characteristics of the reservoir fluid flowing through the flowline 325 and output optical spectra and/or other data representative of the detected optical characteristics. The optical characteristics may include OD of the reservoir fluid at each of the detected wavelengths and/or wavelength ranges. The OD is a logarithmic measurement relating the intensity of light emitted from the light source 315 to the intensity of light detected by the detector 320 at a certain wavelength or range of wavelengths. Each wavelength or wavelength range may correspond to a compositional component of the reservoir fluid. For example, each wavelength or wavelength range may pertain to a corresponding one of CO2, C1, C2, C3, C4, C5, and C6+, although other arrangements are also within the scope of the present disclosure.

The spectrometer 305 may send optical spectra and/or other data representative of the measured optical characteristics to a processor 340 of the control/monitoring system 310. In the context of the present disclosure, the term “processor” refers to any number of processor components. The processor 340 may include a single processor disposed onboard the downhole tool. In other implementations, at least a portion of the processor 340 (e.g., multiple processors collectively operating as the processor 340) may be located within the wellsite surface equipment 160. The processor 340 may also or instead be or include one or more processors located within the downhole tool and connected to one or more processors located in drilling and/or other equipment disposed at the wellsite surface. Moreover, various combinations of processors may be considered part of the processor 340 in the following discussion. Similar terminology is applied with respect to the control/monitoring system 310 as well as a memory 345 of the control/monitoring system 310, meaning that the control/monitoring system 310 may include various processors communicatively coupled to each other and/or various memories at various locations.

The control/monitoring system 310 may estimate or otherwise determine one or more parameters of the reservoir fluid based on the OD data received from the spectrometer 305, a density sensor, a pressure sensor, a temperature sensor, and/or other sensors, and may utilize the parameters to determine density, OD, GOR, mass fractions of compositional components, and/or other properties of the reservoir fluid. To make these and other determinations, the processor 340 may execute instructions stored in the memory 345.

The processor 340 may be communicatively coupled with one or more operator interfaces 350 and/or control devices 355. The operator interface 350 may include logs of predicted and/or measured reservoir fluid properties that are accessible to an operator. The control device 355 may include one or more devices and/or portions thereof that receive control signals for operation based on the estimated properties of the reservoir fluid. Such control devices 355 may implement changes in depth of the downhole tool within the wellbore 111, adjustments to the pumping pressure of the pump(s) of the downhole tool, and/or other control functions, perhaps based on obtained, calculated, and/or estimated reservoir fluid properties.

FIG. 8 is a schematic view of another example wellsite system that can be employed onshore and/or offshore, perhaps including at the same wellsite as depicted in one or more of the figures described above, where the wellbore 111 may have been formed in the one or more subsurface reservoirs by rotary and/or directional drilling. FIG. 8 represents another example environment in which one or more aspects described above may be implemented, such as to perform one or more aspects of the method (50) depicted in FIG. 4 and/or other methods and/or processes within the scope of the present disclosure.

As depicted in FIG. 8, a conveyance means 412 suspended within the wellbore 111 may comprise or be connected to a bottom hole assembly (BHA) 400, which may have a drill bit 405 at its lower end. The conveyance means 412 may comprise drill pipe, wired drill pipe (WDP), tough logging conditions (TLC) pipe, coiled tubing, and/or other means of conveying the BHA 400 within the wellbore 111.

The surface system at the wellsite may comprise a platform, rig, derrick, and/or other wellsite structure 410 positioned over the wellbore 111. The wellsite structure 410 may be substantially similar or identical to the wellsite structure shown in one or more of the figures described above. The wellsite structure 410 may include a rotary table 416, a kelly 417, a hook 418, and/or a rotary swivel 419. The conveyance means 412 may be rotated by the rotary table 416, energized by means not shown, which may engage the kelly 417 at the upper end of the conveyance means 412. The conveyance means 412 may be suspended from the hook 418, which may be attached to a traveling block (not shown), and through the kelly 417 and the rotary swivel 419, which permits rotation of the conveyance means 412 relative to the hook 418. Additionally, or instead, a top drive system may be used.

The surface system may also include drilling fluid 426, which is commonly referred to in the industry as mud, stored in a pit 427 formed at the wellsite. A pump 429 may deliver the drilling fluid 426 to the interior of the conveyance means 412 via a port (not shown) in the swivel 419, causing the drilling fluid to flow downwardly through the conveyance means 412, as indicated by directional arrow 408. The drilling fluid 426 may exit the conveyance means 412 via ports in the drill bit 405, and then circulate upwardly through the annulus region between the outside of the conveyance means 412 and the wall of the wellbore 111, as indicated by directional arrows 409. The drilling fluid 426 may be used to lubricate the drill bit 405 and/or carry cuttings up to the surface as it is returned to the pit 427 for recirculation. Although not pictured, one or more other circulation implementations are also within the scope of the present disclosure, such as a reverse circulation implementation in which the drilling fluid 426 is pumped down the annulus region (i.e., opposite to directional arrows 409) to return to the surface within the interior of the conveyance means 412 (i.e., opposite to directional arrow 408).

The BHA 400 may include various numbers and/or types of downhole tools, schematically depicted in FIG. 8 as downhole tools 420, 430, and 450. One or more of the downhole tools 420, 430, and 450 may be or comprise an acoustic tool, a density tool, a directional drilling tool, a drilling tool, an EM tool, a formation evaluation tool, a gravity tool, a logging while drilling (LWD) tool, a magnetic resonance tool, a measurement while drilling (MWD) tool, a monitoring tool, a neutron tool, a nuclear tool, a photoelectric factor tool, a porosity tool, a reservoir characterization tool, a resistivity tool, a sampling tool, a seismic tool, a surveying tool, a telemetry tool, and/or a tough logging condition (TLC) tool, although other downhole tools are also within the scope of the present disclosure. One or more of the downhole tools 420, 430, and 450 may be utilized to perform at least a portion of a method according to one or more aspects of the present disclosure, perhaps including at least a portion of an implementation of the method (50) shown in FIG. 4.

The downhole tools 420, 430, and/or 450 may be housed in a special type of drill collar, as it is known in the art, and may include capabilities for measuring, processing, and/or storing information, as well as for communicating with the other downhole tools 420, 430, and/or 450, and/or directly with a logging and control system and/or other surface equipment 460. Such communication may utilize one or more conventional and/or future-developed two-way telemetry systems, such as may be or comprise a mud-pulse telemetry system, a wired drill pipe telemetry system, an electromagnetic telemetry system, and/or an acoustic telemetry system, among others within the scope of the present disclosure. One or more of the downhole tools 420, 430, and/or 450 may also comprise an apparatus (not shown) for generating electrical power for use by the BHA 400. Example devices to generate electrical power include, but are not limited to, a battery system and a mud turbine generator powered by the flow of the drilling fluid.

FIG. 9 is a schematic view of at least a portion of apparatus according to one or more aspects of the present disclosure. The apparatus is or comprises a processing system 500 that may execute example machine-readable instructions to implement at least a portion of one or more of the methods and/or processes described herein, and/or to implement a portion of one or more of the example downhole tools, surface equipment, and/or remote equipment described herein, including to perform one or more aspects of the method (50) depicted in FIG. 4 and/or other methods and/or processes within the scope of the present disclosure. For example, the downhole tools 105, 200, 420, 430, and/or 450 shown in FIGS. 5, 6, and 8 may each comprise a corresponding instance of the processing system 500 or a subset of the components of the processing system 500 described below. The remote equipment 101 shown in FIG. 5, the surface equipment 160 shown in FIGS. 5 and 6, the control/monitoring system 310 shown in FIG. 7, and/or the surface equipment 460 shown in FIG. 8 may also each be or comprise a corresponding instance of the processing system 500 or a subset of the components thereof. The processing system 500 may be or comprise, for example, one or more processors, controllers, special-purpose computing devices, servers, personal computers, personal digital assistant (“PDA”) devices, smartphones, internet appliances, and/or other types of computing devices.

The processing system 500 may comprise a processor 512 such as, for example, a general-purpose programmable processor. The processor 512 may comprise a local memory 514, and may execute coded instructions 532 present in the local memory 514 and/or another memory device. The processor 512 may execute, among other things, machine-readable instructions or programs to implement the methods and/or processes described herein. The programs stored in the local memory 514 may include program instructions or computer program code that, when executed by an associated processor, enable the downhole tools 105, 200, 420, 430, and/or 450, the remote equipment 101, the surface equipment 160, the control/monitoring system 310, and/or the surface equipment 460 to perform tasks as described herein. The processor 512 may be, comprise, or be implemented by one or a plurality of processors of various types suitable to the local application environment, and may include one or more of general-purpose computers, special purpose computers, microprocessors, digital signal processors (“DSPs”), field-programmable gate arrays (“FPGAs”), application-specific integrated circuits (“ASICs”), and processors based on a multi-core processor architecture, as non-limiting examples. Of course, other processors from other families are also appropriate.

The processor 512 may be in communication with a main memory, such as may include a volatile memory 518 and a non-volatile memory 520, perhaps via a bus 522 and/or other communication means. The volatile memory 518 may be, comprise, or be implemented by random access memory (RAM), static random access memory (SRAM), synchronous dynamic random access memory (SDRAM), dynamic random access memory (DRAM), RAMBUS dynamic random access memory (RDRAM) and/or other types of random access memory devices. The non-volatile memory 520 may be, comprise, or be implemented by read only memory, flash memory, and/or other types of memory devices. One or more memory controllers (not shown) may control access to the volatile memory 518 and/or the non-volatile memory 520.

The processing system 500 may also comprise an interface circuit 524. The interface circuit 524 may be, comprise, or be implemented by various types of standard interfaces, such as an Ethernet interface, a universal serial bus (USB), a third generation input/output (3GIO) interface, a wireless interface, and/or a cellular interface, among others. The interface circuit 524 may also comprise a graphics driver card. The interface circuit 524 may also comprise a communication device, such as a modem or network interface card, to facilitate exchange of data with external computing devices via a network (e.g., Ethernet connection, digital subscriber line (“DSL”), telephone line, coaxial cable, cellular telephone system, satellite, etc.).

One or more input devices 526 may be connected to the interface circuit 524. The input device(s) 526 may permit a user to enter data and commands into the processor 512. The input device(s) 526 may be, comprise, or be implemented by, for example, a keyboard, a mouse, a touchscreen, a track-pad, a trackball, an isopoint, and/or a voice recognition system, among others.

One or more output devices 528 may also be connected to the interface circuit 524. The output devices 528 may be, comprise, or be implemented by, for example, display devices (e.g., a liquid crystal display or cathode ray tube display (CRT), among others), printers, and/or speakers, among others.

The processing system 500 may also comprise one or more mass storage devices 530 for storing machine-readable instructions and data. Examples of such mass storage devices 530 include floppy disk drives, hard drive disks, compact disk (CD) drives, and digital versatile disk (DVD) drives, among others. The coded instructions 532 may be stored in the mass storage device 530, the volatile memory 518, the non-volatile memory 520, the local memory 514, and/or on a removable storage medium 534, such as a CD or DVD. Thus, the modules and/or other components of the processing system 500 may be implemented in accordance with hardware (embodied in one or more chips including an integrated circuit such as an application specific integrated circuit), or may be implemented as software or firmware for execution by a processor. In particular, in the case of firmware or software, the implementation may be provided as a computer program product including a computer readable medium or storage structure embodying computer program code (i.e., software or firmware) thereon for execution by the processor.

FIG. 10 is a flow-chart diagram of at least a portion of a method (600) according to one or more aspects of the present disclosure. The method (600) may be performed via operation of at least one of a downhole tool, surface equipment, remote equipment, and/or a combination thereof, perhaps including one or more instances of at least a portion of the processing system 500 shown in FIG. 9. The downhole tool may be the downhole tool 105 shown in FIG. 5, the downhole tool 200 shown in FIG. 6, one or more of the downhole tools 420, 430 and 450 shown in FIG. 8, and/or other downhole tools described herein or otherwise within the scope of the present disclosure. The surface equipment may be the surface equipment 160 shown in FIGS. 5 and 6, the surface equipment 460 shown in FIG. 8, and/or other surface equipment within the scope of the present disclosure. The remote equipment may be the remote equipment 101 shown in FIG. 5 and/or other remote equipment within the scope of the present disclosure.

The method (600) includes obtaining (605) data associated with at least one of a subterranean hydrocarbon reservoir, fluid contained within the reservoir, fluid obtained from the reservoir by the downhole tool, and/or a combination thereof. For example, the obtained (605) data may include the example data listed above in Table 1.

Molecular sizes of each of a plurality of compositional components of the fluid contained within the reservoir are then estimated (610) based on the obtained (605) data. For example, the molecular sizes may be estimated (610) utilizing ACGA, and the obtained (605) data utilized with the ACGA to estimate (610) the molecular sizes may include: density, optical density, a bulk fluid solubility parameter, and molar volume of the fluid obtained from the reservoir at each of a plurality of depths beneath the wellsite surface; molar volume, partial density, and a compositional component solubility parameter of at least one of the plurality of compositional components at each of the plurality of depths; and temperature of the reservoir. Estimating (610) the molecular sizes utilizing ACGA may utilize Equation (54) set forth above.

A model of a density inversion of the fluid contained within the reservoir is then generated (615) based on the obtained (605) data and the estimated (610) molecular sizes. The generated (615) density inversion model may assume the fluid contained in the reservoir consists of gas transported deeper into the reservoir via diffusion, maltenes transported deeper into the reservoir via diffusion, and asphaltenes transported deeper into the reservoir via diffusion and expulsion by transported methane. For example, the generated (615) density inversion model may be described by Equations (49)-(53) set forth above, utilizing ones of Equations (1)-(48), or using ones of Equations (1)-(48) directly, among other examples within the scope of the present disclosure.

The method (600) may also comprise utilizing the generated (615) density inversion model to estimate (620) an amount of the density inversion, a depth of the density inversion beneath the wellsite and/or relative to the top of the reservoir, and/or a period of time that has elapsed since the density inversion began to form within the reservoir. The method (600) may also comprise generating (625) a model of a gravity-induced current of the density inversion based on the obtained (605) data and the estimated (620) density inversion amount, depth, and elapsed time. For example, the generated (625) gravity current model may be described by Equation (67) set forth below.

$\begin{matrix} {U_{c} = \frac{\Delta\;\rho\;{gk}}{\phi\mu}} & (67) \end{matrix}$ where:

-   U_(c) is the characteristic gravity current velocity; -   Δρ is the amount of density inversion; -   g is gravitational acceleration; -   k is estimated permeability of the reservoir; -   φ is estimated porosity of the reservoir; and -   μ is measured viscosity of fluid obtained from the reservoir.

However, the generated (625) gravity current model may instead be described by Equation (55) set forth above to account for propagation of the density inversion along a sloped lower boundary of the reservoir.

FIG. 11 is a flow-chart diagram of at least a portion of a method (700) according to one or more aspects of the present disclosure. The method (700) may be performed via operation of at least one of a downhole tool, surface equipment, remote equipment, and/or a combination thereof, perhaps including one or more instances of at least a portion of the processing system 500 shown in FIG. 9. The downhole tool may be the downhole tool 105 shown in FIG. 5, the downhole tool 200 shown in FIG. 6, one or more of the downhole tools 420, 430 and 450 shown in FIG. 8, and/or other downhole tools described herein or otherwise within the scope of the present disclosure. The surface equipment may be the surface equipment 160 shown in FIGS. 5 and 6, the surface equipment 460 shown in FIG. 8, and/or other surface equipment within the scope of the present disclosure. The remote equipment may be the remote equipment 101 shown in FIG. 5 and/or other remote equipment within the scope of the present disclosure. At least a portion of the method (700) may also be performed in conjunction with at least a portion of the method (600) shown in FIG. 10.

The method (700) includes obtaining (705) data associated with at least one of a subterranean hydrocarbon reservoir, fluid contained within the reservoir, fluid obtained from the reservoir by the downhole tool, and/or a combination thereof. For example, the obtained (705) data may include the example data listed above in Table 1.

The method (700) also includes identifying (710) at least one condition indicating that a density inversion exists in the fluid contained in the reservoir. For example, the fluid obtained from the reservoir may be via operation of the downhole tool at each of a plurality of depths within the reservoir beneath the wellsite, and the obtained (705) data may include GOR data related to the fluid obtained at each of the plurality of depths. In such implementations, the at least one identified (710) condition may be a disequilibrium gradient of GOR relative to depth within the reservoir beneath the wellsite, as determined based on the GOR data at each of the plurality of depths.

In another example, the fluid obtained from the reservoir may again be via operation of the downhole tool at each of a plurality of depths within the reservoir beneath the wellsite, but the obtained (705) data may include composition data related to the fluid obtained at each of the plurality of depths, such as a mass concentration of asphaltenes at each of the plurality of depths. In such implementations, the at least one identified (710) condition may be an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths.

In another example, the fluid obtained from the reservoir may again be via operation of the downhole tool at each of a plurality of depths within the reservoir beneath the wellsite, and the obtained (705) data may include both GOR and composition data related to the fluid obtained at each of the plurality of depths. Accordingly, the at least one identified (710) condition includes at least one of: a first condition in which there exists a disequilibrium gradient of GOR relative to depth within the reservoir beneath the wellsite, as determined based on the GOR data at each of the plurality of depths; and a second condition in which there exists an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths. In such implementations, the method (700) may also comprise identifying (715) a current status of a gravity-induced current of asphaltenes towards the lower base of the reservoir as a result of the formation of the density inversion. For example, if both the first and second conditions exist, then identifying (715) the current status of the gravity current may include identifying that the density inversion is presently causing an ongoing, gravity-induced current of asphaltenes towards the lower base of the reservoir. However, if the second condition exists, but the first condition does not exist, then identifying (715) the current status of the gravity current may include identifying that the density inversion previously caused a now-substantially-completed, gravity-induced current of asphaltenes towards the lower base of the reservoir.

In view of the entirety of the present disclosure, including the claims and the figures, a person having ordinary skill in the art will readily appreciate that the present disclosure introduces a system comprising: a downhole tool operable within a wellbore extending from a wellsite surface into a subterranean hydrocarbon reservoir, wherein the downhole tool is operable to obtain fluid from the reservoir and data associated with the obtained fluid; and an apparatus operable to identify at least one condition indicating that a density inversion exists in the fluid contained in the reservoir based on the obtained data.

The downhole tool may be operable to obtain the fluid from the reservoir at each of a plurality of depths within the reservoir, the obtained data may include GOR data related to the fluid obtained at each of the plurality of depths, and the at least one condition may be a disequilibrium gradient of GOR relative to depth within the reservoir, as determined based on the GOR data at each of the plurality of depths.

The downhole tool may be operable to obtain the fluid from the reservoir at each of a plurality of depths within the reservoir, the obtained data may include composition data related to the fluid obtained at each of the plurality of depths, and the at least one condition may be an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths.

The downhole tool may be operable to obtain the fluid from the reservoir at each of a plurality of depths within the reservoir, the obtained data may include GOR and composition data related to the fluid obtained at each of the plurality of depths, and the at least one condition may include at least one of: a first condition in which there exists a disequilibrium gradient of GOR relative to depth within the reservoir, as determined based on the GOR data at each of the plurality of depths; and a second condition in which there exists an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths. In such implementations, among others within the scope of the present disclosure, the apparatus may be further operable to (a) identify, based on the existence of both of the first and second conditions, that the density inversion is causing an ongoing, gravity-induced current of asphaltenes towards the lower base of the reservoir; (b) to identify, based on the existence of the second condition and the absence of the first condition, that the density inversion previously caused a now-substantially-completed, gravity-induced current of asphaltenes towards the lower base of the reservoir; or (c) identify, based on the existence of both of the first and second conditions, that the density inversion is causing an ongoing, gravity-induced current of asphaltenes towards the lower base of the reservoir; and identify, based on the existence of the second condition and the absence of the first condition, that the density inversion previously caused a now-substantially-completed, gravity-induced current of asphaltenes towards the lower base of the reservoir.

The apparatus may be further operable to: estimate molecular sizes of each of a plurality of compositional components of the fluid contained within the reservoir; and generate a model of the density inversion of the fluid contained within the reservoir based on the obtained data and the estimated molecular sizes. In such implementations, among others within the scope of the present disclosure, the apparatus may be further operable to, based on the generated density inversion model, estimate an amount of the density inversion, a depth of the density inversion, and a period of time elapsed since the density inversion began to form within the reservoir. The apparatus may be further operable to generate a model of a gravity-induced current of the density inversion based on the obtained data and the estimated density inversion amount, depth, and elapsed time. The apparatus may be operable to estimate the molecular sizes utilizing an asphaltene concentration gradient analysis (ACGA) and the obtained data. For example, the obtained data utilized with the ACGA to estimate the molecular sizes may include: each of density, optical density, a bulk fluid solubility parameter, and molar volume of the fluid obtained from the reservoir at each of a plurality of depths within the reservoir; each of molar volume, partial density, and a compositional component solubility parameter of at least one of the plurality of compositional components at each of the plurality of depths; and temperature of the reservoir. The generated density inversion model may assume the fluid contained in the reservoir consists of: gas transported deeper into the reservoir via diffusion; maltenes transported deeper into the reservoir via diffusion; and asphaltenes transported deeper into the reservoir via diffusion and expulsion by transported gas.

Inputs for the density inversion model may include: pressure of the reservoir; temperature of the reservoir; density of the fluid contained in the reservoir; viscosity of the fluid contained in the reservoir; composition of the fluid contained in the reservoir; depth within the reservoir; time; the estimated molecular sizes; a diffusion coefficient DGA between gas and asphaltene components of the fluid contained in the reservoir; a diffusion coefficient DGM between gas and maltene components of the fluid contained in the reservoir; and a diffusion coefficient DAM between asphaltene and maltene components of the fluid contained in the reservoir. The downhole tool may be operable to obtain ones of the density inversion model inputs, including: the pressure of the reservoir; the temperature of the reservoir; the density of the fluid contained in the reservoir; the viscosity of the fluid contained in the reservoir; and the composition of the fluid contained in the reservoir.

The present disclosure also introduces a system comprising: (a) a downhole tool operable to obtain: (i) fluid from a subterranean reservoir at each of a plurality of depths within the reservoir; and (ii) data associated with the obtained fluid, including gas-oil ratio (GOR) and composition data related to the fluid obtained at each of the plurality of depths; and (b) an apparatus operable to identify at least one condition indicating that a density inversion exists in the fluid contained in the reservoir based on the obtained data, wherein the at least one condition includes at least one of: (i) a first condition in which there exists a disequilibrium gradient of GOR relative to depth within the reservoir, as determined based on the GOR data at each of the plurality of depths; and (ii) a second condition in which there exists an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths.

The apparatus may be further operable to: identify, based on the existence of both of the first and second conditions, that the density inversion is causing an ongoing, gravity-induced current of asphaltenes towards the lower base of the reservoir; and identify, based on the existence of the second condition and the absence of the first condition, that the density inversion previously caused a now-substantially-completed, gravity-induced current of asphaltenes towards the lower base of the reservoir.

The apparatus may be further operable to: estimate molecular sizes of each of a plurality of compositional components of the fluid contained within the reservoir; generate a model of the density inversion of the fluid contained within the reservoir based on the obtained data and the estimated molecular sizes, wherein the generated density inversion model may assume the fluid contained in the reservoir consists of gas, maltenes, and asphaltenes; and based on the generated density inversion model, estimate an amount of the density inversion, a depth of the density inversion, and a period of time elapsed since the density inversion began to form within the reservoir. The apparatus may be further operable to generate a model of a gravity-induced current of the density inversion based on the obtained data and the estimated density inversion amount, depth, and elapsed time.

The present disclosure also introduces a method comprising: (1) operating a downhole tool to obtain: (a) fluid from a subterranean reservoir at each of a plurality of depths within the reservoir; and (b) data associated with the obtained fluid, including gas-oil ratio (GOR) and composition data related to the fluid obtained at each of the plurality of depths; and (2) operating an apparatus comprising a processor and a memory including computer program code, including to: (a) identify at least one condition indicating that a density inversion exists in the fluid contained in the reservoir based on the obtained data, wherein the at least one condition includes at least one of: (i) first condition in which there exists a disequilibrium gradient of GOR relative to depth within the reservoir, as determined based on the GOR data at each of the plurality of depths; and (ii) a second condition in which there exists an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths; (b) identify, based on the existence of both of the first and second conditions, whether the density inversion is causing an ongoing, gravity-induced current of asphaltenes towards the lower base of the reservoir; (c) identify, based on the existence of the second condition and the absence of the first condition, whether the density inversion previously caused a now-substantially-completed, gravity-induced current of asphaltenes towards the lower base of the reservoir; (d) estimate molecular sizes of each of a plurality of compositional components of the fluid contained within the reservoir utilizing an asphaltene concentration gradient analysis (ACGA) and the obtained data; (e) generate a model of the density inversion of the fluid contained within the reservoir based on the obtained data and the estimated molecular sizes, wherein the generated density inversion model assumes the fluid contained in the reservoir consists of gas, maltenes, and asphaltenes; (f) based on the generated density inversion model, estimate an amount of the density inversion, a depth of the density inversion, and a period of time elapsed since the density inversion began to form within the reservoir; and (g) generate a model of a gravity-induced current of the density inversion based on the obtained data and the estimated density inversion amount, depth, and elapsed time.

The foregoing outlines features of several embodiments so that a person having ordinary skill in the art may better understand the aspects of the present disclosure. A person having ordinary skill in the art should appreciate that they may readily use the present disclosure as a basis for designing or modifying other processes and structures for carrying out the same purposes and/or achieving the same advantages of the embodiments introduced herein. A person having ordinary skill in the art should also realize that such equivalent constructions do not depart from the scope of the present disclosure, and that they may make various changes, substitutions and alterations herein without departing from the spirit and scope of the present disclosure.

The Abstract at the end of this disclosure is provided to comply with 37 C.F.R. § 1.72(b) to allow the reader to quickly ascertain the nature of the technical disclosure. It is submitted with the understanding that it will not be used to interpret or limit the scope or meaning of the claims. 

What is claimed is:
 1. A system comprising: a downhole tool operable within a wellbore extending from a wellsite surface into a subterranean hydrocarbon reservoir, wherein the downhole tool is operable to obtain fluid from the reservoir and data associated with the obtained fluid; and an apparatus operable to identify at least one condition indicating that a density inversion exists in the fluid contained in the reservoir based on the obtained data, wherein the apparatus is operable to: estimate molecular sizes of each of a plurality of compositional components of the fluid contained within the reservoir; and generate a model of the density inversion of the fluid contained within the reservoir based at least in part on the obtained data and the estimated molecular sizes, wherein inputs for the density inversion model include one or more diffusion coefficients related to gas, maltene, asphaltene, or any combination thereof.
 2. The system of claim 1 wherein: the downhole tool is operable to obtain the fluid from the reservoir at each of a plurality of depths within the reservoir; the obtained data includes gas-oil ratio (GOR) data related to the fluid obtained at each of the plurality of depths; and the at least one condition is a disequilibrium gradient of GOR relative to depth within the reservoir, as determined based on the GOR data at each of the plurality of depths.
 3. The system of claim 1 wherein: the downhole tool is operable to obtain the fluid from the reservoir at each of a plurality of depths within the reservoir; the obtained data includes composition data related to the fluid obtained at each of the plurality of depths; and the at least one condition is an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths.
 4. The system of claim 1 wherein: the downhole tool is operable to obtain the fluid from the reservoir at each of a plurality of depths within the reservoir; the obtained data includes gas-oil ratio (GOR) and composition data related to the fluid obtained at each of the plurality of depths; and the at least one condition includes at least one of: a first condition in which there exists a disequilibrium gradient of GOR relative to depth within the reservoir, as determined based on the GOR data at each of the plurality of depths; and a second condition in which there exists an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths.
 5. The system of claim 4 wherein the apparatus is further operable to identify, based on the existence of both of the first and second conditions, that the density inversion is causing an ongoing, gravity-induced current of asphaltenes towards the lower base of the reservoir.
 6. The system of claim 4 wherein the apparatus is further operable to identify, based on the existence of the second condition and the absence of the first condition, that the density inversion previously caused a now-substantially-completed, gravity-induced current of asphaltenes towards the lower base of the reservoir.
 7. The system of claim 4 wherein the apparatus is further operable to: identify, based on the existence of both of the first and second conditions, that the density inversion is causing an ongoing, gravity-induced current of asphaltenes towards the lower base of the reservoir; and identify, based on the existence of the second condition and the absence of the first condition, that the density inversion previously caused a now-substantially-completed, gravity-induced current of asphaltenes towards the lower base of the reservoir.
 8. The system of claim 1 wherein the apparatus is further operable to estimate, based on the generated density inversion model, an amount of the density inversion, a depth of the density inversion, and a period of time elapsed since the density inversion began to form within the reservoir.
 9. The system of claim 8 wherein the apparatus is further operable to generate a model of a gravity-induced current of the density inversion based on the obtained data and the estimated density inversion amount, depth, and elapsed time.
 10. The system of claim 1 wherein the apparatus is operable to estimate the molecular sizes utilizing an asphaltene concentration gradient analysis (ACGA) and the obtained data.
 11. The system of claim 10 wherein the obtained data utilized with the ACGA to estimate the molecular sizes includes: each of density, optical density, a bulk fluid solubility parameter, and molar volume of the fluid obtained from the reservoir at each of a plurality of depths within the reservoir; each of molar volume, partial density, and a compositional component solubility parameter of at least one of the plurality of compositional components at each of the plurality of depths; and temperature of the reservoir.
 12. The system of claim 1 wherein the generated density inversion model assumes the fluid contained in the reservoir consists of: gas transported deeper into the reservoir via diffusion; maltenes transported deeper into the reservoir via diffusion; and asphaltenes transported deeper into the reservoir via diffusion and expulsion by transported gas.
 13. The system of claim 1 wherein the inputs for the density inversion model include: pressure of the reservoir; temperature of the reservoir; density of the fluid contained in the reservoir; viscosity of the fluid contained in the reservoir; composition of the fluid contained in the reservoir; depth within the reservoir; time; the estimated molecular sizes; and wherein the one or more diffusion coefficients include: a diffusion coefficient D_(GA) between gas and asphaltene components of the fluid contained in the reservoir; a diffusion coefficient D_(GM) between gas and maltene components of the fluid contained in the reservoir; and a diffusion coefficient D_(AM) between asphaltene and maltene components of the fluid contained in the reservoir.
 14. The system of claim 13 wherein the downhole tool is operable to obtain ones of the density inversion model inputs, including: the pressure of the reservoir; the temperature of the reservoir; the density of the fluid contained in the reservoir; the viscosity of the fluid contained in the reservoir; and the composition of the fluid contained in the reservoir.
 15. A system comprising: a downhole tool operable to obtain: fluid from a subterranean reservoir at each of a plurality of depths within the reservoir; and data associated with the obtained fluid, including gas-oil ratio (GOR) and composition data related to the fluid obtained at each of the plurality of depths; and an apparatus operable to identify at least one condition indicating that a density inversion exists in the fluid contained in the reservoir based on the obtained data, wherein the at least one condition includes at least one of: a first condition in which there exists a disequilibrium gradient of GOR relative to depth within the reservoir, as determined based on the GOR data at each of the plurality of depths; and a second condition in which there exists an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths; wherein the apparatus is operable to: estimate molecular sizes of each of a plurality of compositional components of the fluid contained within the reservoir; and generate a model of the density inversion of the fluid contained within the reservoir based at least in part on the obtained data and the estimated molecular sizes, wherein the generated density inversion model assumes the fluid contained in the reservoir comprises gas, maltenes, and asphaltenes, and wherein inputs for the model of the density inversion comprise one or more diffusion coefficients related to gas, maltene, and asphaltene.
 16. The system of claim 15 wherein the apparatus is further operable to: identify, based on the existence of both of the first and second conditions, that the density inversion is causing an ongoing, gravity-induced current of asphaltenes towards the lower base of the reservoir; and identify, based on the existence of the second condition and the absence of the first condition, that the density inversion previously caused a now-substantially-completed, gravity-induced current of asphaltenes towards the lower base of the reservoir.
 17. The system of claim 15 wherein the apparatus is further operable to: based on the generated density inversion model, estimate an amount of the density inversion, a depth of the density inversion, and a period of time elapsed since the density inversion began to form within the reservoir.
 18. The system of claim 17 wherein the apparatus is further operable to generate a model of a gravity-induced current of the density inversion based on the obtained data and the estimated density inversion amount, depth, and elapsed time.
 19. A method comprising: operating a downhole tool to obtain: fluid from a subterranean reservoir at each of a plurality of depths within the reservoir; and data associated with the obtained fluid, including gas-oil ratio (GOR) and composition data related to the fluid obtained at each of the plurality of depths; and operating an apparatus comprising a processor and a memory including computer program code, including to: identify at least one condition indicating that a density inversion exists in the fluid contained in the reservoir based on the obtained data, wherein the at least one condition includes at least one of: a first condition in which there exists a disequilibrium gradient of GOR relative to depth within the reservoir, as determined based on the GOR data at each of the plurality of depths; and a second condition in which there exists an accumulation of asphaltenes proximate a lower base of the reservoir, as determined based on the composition data at each of the plurality of depths; identify, based on the existence of both of the first and second conditions, whether the density inversion is causing an ongoing, gravity-induced current of asphaltenes towards the lower base of the reservoir; identify, based on the existence of the second condition and the absence of the first condition, whether the density inversion previously caused a now-substantially-completed, gravity-induced current of asphaltenes towards the lower base of the reservoir; estimate molecular sizes of each of a plurality of compositional components of the fluid contained within the reservoir utilizing an asphaltene concentration gradient analysis (ACGA) and the obtained data; generate a model of the density inversion of the fluid contained within the reservoir based at least in part on the obtained data and the estimated molecular sizes, wherein the generated density inversion model assumes the fluid contained in the reservoir comprises gas, maltenes, and asphaltenes and wherein inputs for the model of the density inversion comprise one or more diffusion coefficients related to gas, maltene, asphaltene, or any combination thereof; based on the generated density inversion model, estimate an amount of the density inversion, a depth of the density inversion, and a period of time elapsed since the density inversion began to form within the reservoir; and generate a model of a gravity-induced current of the density inversion based on the obtained data and the estimated density inversion amount, depth, and elapsed time.
 20. The system of claim 15, wherein the generated density inversion model assumes the fluid contained in the reservoir consists of gas, maltenes, and asphaltenes. 